Peristaltic modes of single vortex in the Abelian Higgs model 
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Using the Abelian Higgs model, we study the radial excitations of single vortex and their propa- 
gation modes along the vortex line. We call such beyond-stringy modes peristaltic modes of single 
vortex. With the profile of the static vortex, we derive the vortex-induced potential, i.e., single- 
particle potential for the Higgs and the photon field fluctuations around the static vortex, and 
investigate the coherently propagating fluctuations which corresponds to the vibration of the vor- 
tex. We derive, analyze and numerically solve the field equations of the Higgs and the photon field 
fluctuations around the static vortex with various Ginzburg-Landau parameter k and topological 
charge n. Around the BPS value or critical coupling fi:^ = l/2, there appears a signiflcant correlation 
between the Higgs and the photon fleld fluctuations mediated by the static vortex. As a result, 
for K^ = l/2, we find the characteristic new-type discrete pole of the peristaltic mode corresponding 
to the quasi-bound-state of coherently fiuctuating fields and the static vortex. We investigate its 
excitation energy, correlation energy of coherent fluctuations, spatial distributions, and the resulting 
magnetic flux behavior in detail. Our investigation covers not only usual Type-H vortices with n=l 
but also Type-I and Type-H vortices with n £ Z for the application to various general systems 
where the vortex-like objects behave as the essential degrees of freedom. 

PACS numbers: 11.27.-|-d, ll.15.Kc, 74.25.-q, 12.38-t 



> 

^ • 
o ■ 
(N : 
o . 
I> ■ 
o 

Oh. 
I 



X 



I. INTRODUCTION 

There have been wide interests in the study of the topo- 
logical objects such as the vortices monopoles 0,11], 
and instantons [1] because of their importance for under- 
standing some crucial aspects of the nonlinear field the- 
ories 5, 6]. They are not included in the original action 
manifestly, but appear as the localized finite energy solu- 
tions of the nonlinear field equations. These solutions can 
be classified with the topological number which comes 
from the topological property of the system. Since the 
topological objects have spatially localized, finite energy 
configurations and some topologically conserved charge, 
it is fascinating to regard them as dynamical degrees of 
freedom in the system. 

One of such objects is the Abrikosov vortex 1] in 
Ginzburg-Landau (GL) theory, which appears as the 
magnetic flux squeezed by the Cooper-pair condensate. 
When the external magnetic field applied to the super- 
conductor exceeds some critical strength, the Cooper- 
pair is dissociated and the node of the Cooper-pair forms 
a vortex line. The single- valued property of the Cooper- 
pair around vortex lines leads the quantization of the 
total magnetic flux, as e$ = 27rrt {n E Z), where e, <i>,n 
represent the effective gauge coupling of the Cooper-pair, 
the total magnetic flux and the topological number, re- 
spectively. The topological number n characterizes the 
topology of the system and we can classify vortices in 
terms of n. 

There exist many applications employing the concept 
of the Abrikosov vortex and its relativistic version, the 
Nielsen-Olesen vortex j7|. One example is the squeezed 
color-electric flux between the quark and anti- qua rk in 
QCD. Nambu [|, 't Hooft and Mandelstam [13] dis- 
cussed that there exist the vortices which appear as the 



color-electric flux tubes squeezed by the condensation of 
magnetic monopoles. As a result of the color-electric 
flux squeezing, we can only observe colorless mesons and 
baryons, which are color singlet bound states of the col- 
ored quarks with the bond of the color-electric flux. This 
scenario is a possible explanation for "color confinement" , 
which is the experimental fact in hadron physics pd| . 
The lattice QCD Monte Carlo calculations dz] for QQ 
[Tst and 3Q [3, [11] potentials strongly support the flux 
tube picture, and show the universality of the string ten- 
sion of the flux-tube, a ~ 0.89 GeV, in both QQ and 3Q 
cases. 

An another example is the cosmic string in astrophysics 
as a seed of the galaxy formation and baryogensis [l^. 
Kibble [13] and Zurek [l^ discussed the distribution of 
the topological defects and their cosmological evolution 
after the quench of the early Universe and the resulting 
phase transition. They argued that the cosmic strings 
as topological defects would produce a characteristic sig- 
nature in the cosmic microwave background. Following 
their scenarios, the phase transition dynamics are closely 
investigated in the numerical simulations using the phe- 
nomenological time-dependent GL (TDGL) theory in 
both cases with global symmetry and local gauge symme- 
try. In particular, in the case of the gauge symmetry, it is 
discussed that the initial thermal magnetic fiuctuations 
just after the quench play important roles on the later 
distribution of the vortices such as the clusters of vor- 
tices [H, [2l| . These investigations are closely correlated 
with the condensed matter physics in a viewpoint to un- 
derstand the phase transition dynamics, and provide the 
interesting subjects to link different fields of research. 

Now we turn to the dynamical aspects of single vortex 
as a "pseudo-particle" . For the treatment of the vortex 
motion, the vortex is usually regarded as one dimensional 
(1-D) thin object like a string, neglecting the extension 
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perpendicular to the vortex line. Such a treatment is ex- 
pected to be sufRcient when the length of the string is 
much larger than the length scale in the radial direction. 
This condition comes from the fact that the excitation 
energy of string as a 1-D object is typically ^ tt/Is {Is 
is the string length). Then, for the vortex with large 
length, the stringy excitation is more important than ra- 
dial one with changing its thickness. Even if the above 
condition is not well-satisfied, we have only to consider 
extensive modes of the magnetic flux when the Type-II 
nature is very strong, i.e., the Ginzburg-Landau parame- 
ter K = i5/^ is very large (Here 5 is the penetration depth 
of the magnetic field and ^ is the coherence length of 
the Higgs field). This is because when Type-II nature 
is very strong, the topological defect of the Higgs field 
along the vortex line is strongly squeezed and its exten- 
sion modes perpendicular to the vortex line cost a large 
energy. Then, in the Type-II limit of the vortex, the main 
excitation is given by the stringy modes, and the subject 
is linking to one of the most fundamental problems in 
the quantum string theory [S^], which gives not only a 
candidate of the grand- unification including gravity but 
also a new method to analyze nonperturbative QCD 23]. 
So far, there are so many interesting studies done on the 
quantum string. For instance, the quantum fluctuations 
of the string lead to a "roughening", i.e., a widening of 
the string 24], and several studies of quantum Niclscn- 
Olesen strings [25l [2^ showed that quantum vibrational 
modes make the vortex world-sheets crumpled and dom- 
inated by "branched-polymers" as a disease of quantum 
strings [25|. 

However, we sometimes encounter the cases where the 
vortex cannot be regarded as a stringy object. We are 
interested in such cases. One example is the Type-II su- 
perconductors in condensed matter physics near k^=1/2, 
i.e., Bogomol'nyi-Prasad-Sommerfeld (BPS) [13, HI] 
value or critical coupling. In such cases, the squeezing 
of the magnetic flux due to Cooper-pair condensation is 
not very strong and both the Cooper-pair and photon 
fluctuations can play important roles. Moreover, their 
interplay may provide interesting phenomena. 

An another interesting example is the quark systems 
such as the meson and baryon which have the color- 
electric flux tube with the total length about I fm, which 
is not very large in comparison with the flux tube ex- 
tension of about 0.4 fm. In addition, accordin g to the 
analysis using the dual-Ginzburg-Landau model [29l [sot , 
Type-II nature also seems to be not so strong, and then 
the radial motion of the vortex can play an important 
role. Studies on such excitation modes are important to 
understand the structure of hadrons. For the gluonic ex- 
citations, the lattice QCD studies for heavy QQ 'sT] and 
3Q [13] system both provide the 1st gluonic excitation en- 
ergy as Ai? ^ 1 GeV in the the total length of the flux, 
-Z^min=0. 5-1.0 fm, as the typical hadronic scale. Small de- 
pendence on Lmin provides the possibility to understand 
the 1st gluonic excitation as not simple stringy vibra- 
tion but the radial excitation with the fixed boundary 



effects on the flux-tube, because the simple stringy mode 
is likely to behave proportionally to l/Lmin- The studies 
on the energy of the radial excitation with comparing it 
to that of lattice results will give a qualitative explana- 
tion for the gluonic excitations in the case of relatively 
small but not too short Lmin- There seem to exist several 
other cases where the radial excitations play important 
roles. 

Putting applications to the above examples in perspec- 
tive, we will discuss the axial-symmetric radial excitation 
modes and their propagation along the vortex line. Here- 
after, we will call this excitation mode peristaltic mode. 
The schematic picture of the peristaltic mode is shown 
in Fig[TJ To discuss the peristaltic modes, we utilize the 
Abelian Higgs action which includes the Higgs and the 
photon fields. In terms of this action, the peristaltic 
modes are characterized by the Higgs and photon field 
fluctuations. In particular, the time-dependent photon 
fluctuations, which are sometimes neglected, play impor- 
tant roles on the corporative behavior with the Higgs field 
fiuctuations. These fluctuations are regarded as the vi- 
bration of the vortex from the viewpoint that the vortex 
is the dynamical degrees of freedom in the system. 

In this work, since we concentrate on the radial excita- 
tion modes, we simply neglect the effects of the boundary 
of the vortex edge, considering the vortex with the infl- 
nite length. The bending of the vortex is also neglected. 
We investigate both Type-I and Type-II vortices with 
not only usual n=l but also |n| > 2 for the future appli- 
cations to the non-equilibrium or multi-vortices systems 
where the giant vortices can appear as excited states. In 
the Type-I case, if the system contains only flnite mag- 
netic flux in contrast to the equilibrium case usually re- 
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FIG. 1: The illustration for the peristaltic mode of the vortex. 
The axial-symmetric oscillation perpendicular to the vortex 
line propagates with the excitation of the Higgs and photon 
fields. The waving solid line in the left side represents the 
surface where the amplitude of the Higgs field equals some 
constant value. The right hand side represents the amplitude 
of the Higgs and the magnetic fields in the cross section per- 
pendicular to the vortex line. Through this process, the con- 
served total magnetic flux repeatedly contracts and relaxes. 
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alized in the laboratory, these flux gather and become a 
|n| > 2 giant vortex because Type-I vortices attract each 
other. On the other hand, in the Type-II case, the gi- 
ant vortex is unstable and splits into the small vortices. 
In spite of this, when is not too large, it is worth 
considering the vibration of the giant vortex because the 
gradient of the vortex-vortex interaction is small for the 
small separation of vortices 33] and we can expect that 
the fission process develops slowly. 

The organization of this paper is as follows. In Sec. II 
we briefly review the treatments of the static profile of 
single vortex. In Sec. Ill we derive the effective action 
and the resulting equation of motion for the fluctuations 
around the vortex. The effects of the static vortex back- 
ground on the fluctuations are emphasized. In Sec. IV we 
show the numerical results for the excitation modes of 
the n=l vortex with 0.1<k^<1.0. Based on the classi- 
fication with K^, we closely discuss their excitation en- 
ergy, correlation energy between the Higgs and photon 
fields, spatial distribution of fluctuations, and the result- 
ing magnetic flux behavior. It will be shown that, in the 
case of ~ 1/2, the interplay of the Higgs and the pho- 
ton fields plays an essential role for the existence of the 
new-type discrete pole. We also discuss the excitation 
modes around \n\ > 2 vortices. In Sec.V we summarize 
this paper with the perspective for the future direction. 



II. ABELIAN HIGGS MODEL AND THE 
STATIC VORTEX 

In this section, we briefly review the treatments of the 
static vortex in terms of Abelian Higgs model [7|. This 
type of action is very popular and has been used to de- 
scribe cosmic string in the astrophysics 16], the color flux 
tube in QCD [33] , and so on. We use this Lorentz invari- 
ant action to investigate the properties of the vortex as 
a "pseudo-particle". 

For the finite temperature cases, one usually uses the 
phenomenological model with replacing the 2nd time- 
derivative in Abelian Higgs model with the 1st derivative 
although the microscopic foundation of this treatment is 
still not conclusive. The noise term is sometimes fur- 
ther added. These replacements change the equation of 
motion from the wave equation to the diffusive equation 
which respects the dissipation effects due to the thermal 
fluctuations. We will discuss the finite temperature case 
elsewhere. 



A. The action and the field equation 

Let us consider the ideal system which contains only 
one vortex with infinite length along the ^-axis. Our pur- 
pose is to investigate the excitation mode of the vortex 
in such a system. We start with the following action with 



natural units {c = h = 1) 

-x{\^\'-vr], (1) 

where ip and A^ represent the Higgs field and photon field 
respectively. Here e, A, v are the effective gauge coupling 
constant, the strength of self-coupling, and the vacuum 
expectation value of the Higgs field, respectively, x'^ de- 
note space-time coordinates and 9'^ = Through this 
paper, we describe the variables in physical unit with bar. 
For later convenience, we rescale the fields and the coor- 
dinates to make them dimensionless quantities, 

iP = iP/v, A^' = A'^/v, x^' = era^. (2) 

We also rewrite the action 

S = e^S. (3) 

The energy relation between the physical and our 
rescaled unit can be obtained with noting the following 
relation, 

S = J dt L = J dt L ^ J dt ^L. (4) 

Here L and L are Lagrangian in physical and rescaled 
unit, respectively, and they are related as Z = Sim- 
ilarly, we can obtain the energy and momentum in the 
physical unit with multiplying v /e to our rescaled energy 
and momenta. 

After these rewriting, the action reads 

S = J d^xS^ 1 {d" - 1^ -^{d'^A'' ~ d^'A^f 

(5) 

Here we have introduced the GL parameter k} — Xj — 
6/£_, which is used to distinguish the Type-I (k < l/\/2 ) 
and Type-II (k > 1/V2 ) superconductors. The critical 
case of K = l/\/2 is called as the BPS case, which lies 
between Type-I and Type-II. As a remarkable fact in 
the BPS case, the photon mass tua and the Higgs mass 
7710 coincide as niA = rn^^ because of mA — \/2ev and 
rrifj, = 2^/\v, i.e., 2k^ = 2A/e^ = m'^/m\, which can be 
derived from the original action ([T]). 

In this work, we will consider the system which is trans- 
lation invariant in the z-direction and contains single vor- 
tex with infinite length along the z-axis. We also impose 
the axial-symmetry around the z-axis on the static so- 
lution. With cyhndrical coordinates {r,6,z), we denote 
the complex field tpix) in the polar-decomposed form 
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where n is the topological number (n=0, ±1, . . . ) related 
to the phase of the Higgs field wave function around the 
z-axis. For single vortex with the topological number 
n S Z, the sign of n is not relevant because of charge 
conjugation symmetry of the Abelian Higgs model, and 
therefore we have only to consider n=0, 1, . . . cases with- 
out loss of generality. Since the Abelian Higgs model has 
the U(l) local gauge symmetry, we have the degrees of 
freedom to choose the gauge. For the treatment of the 
vortex solution, the familiar gauge is the Aq = gauge. 
Instead, we adopt the gauge which eliminates the small 
amplitude of the phase. 



X = 0, 



(7) 



to simplify the equations of fluctuations in the later anal- 
ysis. For n—0, this gauge fixing corresponds to the uni- 
tary gauge. 

From the rotational symmetry and translational invari- 
ance along the z-direction, we adopt the following Ansatz 
for the static solutions 

= (0, -Ag{r) sine, Ag{r) cosO, ), (8) 

and the static electric and magnetic fields are obtained 
as 

E(r)=0, H(r) = (0, 0, -^{rAe) ). (9) 

r ar 



As seen soon later in Eq. p3)l and (|T4)l . our gauge fixing 
and Ansatz for the static fields reproduce the same results 
for the static field solutions as those in usual Aq — 
gauge. 

The total static energy of the system (in the rescaled 
unit) is obtained as E = —S with putting the time deriva- 
tive terms zero, 



-^static = jdzdr 2'KrE[4>,Ag], 



(10) 



where 



+ \{d,AgY + ]^{drAe + 



[n-rAef 



(The full action and the energy functional are given in 
Appendix Under the above circumstances, we con- 
sider the static vortex solutions (j)(r) and Ag{r) which 
minimize the static energy. From the variational princi- 
ple 

^-^static ^-^static ^ /lo\ 

5(p{r) ~ 5Ae{r) ~ ' ^ ' 

and then the static field equations read 

(r^) + ^ :r-^(l) + 2K^(l>{(f>^ -1) =0, (13) 



r dr dr 

d / dAfi \ Ao , , N , o 
-r {r-r-) - ^ + 2 n - rA^ = 0. 
ar ar r 



(14) 



Since we are interested in the vortex solution, we consider 
the boundary condition (B.C.) which is appropriate to 
describe the vortex at the center. For r ^ 0, 



(r) const. X r", Ag{r) const, x r. 



(15) 



and for 



(r) — > 1 -I- const. X r ^e' 
n 



'2Kr 



Ae(r) ^ -+ const. X r . 
r 



(16) 



These boundary conditions are obtained by the asymp- 
totic analysis of the static field equations which satisfy 
the following physical requirements to describe the vor- 
tex at the center: (I) at r = the vacuum expectation 
value of Higgs field takes zero. (II) the static field en- 
ergy U asymptotically approaches to the vacuum energy 
when r — s- cxD. Note that the first condition represents 
the existence of the vortex at r = and the second one 
represents the magnetic fields are completely screened in 
the region far from the center of the vortex. The bound- 
ary conditions (fT5)) and (fTB|) also lead the quantization 
of the total magnetic flux $0i 



r°° f°° d 

$0 = / dr 2TTrH^{r) = / dr 2Tr—{rAe{r)) 

Jo Jo dr 

= 2Tm. 



(17) 



(<I>o corresponds to e$ in the physical unit.) 

We can numerically solve the static field equations (fT3)) 
and (fT4|l with the boundary conditions (fT5|) and (fT6|) and 
obtain the static solutions (/)'^'(r) and Af{r) for the Higgs 
field (l){r) and the photon field Ag (r) . Figl^] shows the 
static profiles of the Higgs field ^i'^'(r), the magnetic field 
H'^^{r) = ■i^(rAg'), and the resulting static energy dis- 
tribution £{r) in the vortex with various n. 



B. The vortex mass and the classification of Type-I 
and Type-II superconductors 

Using the distribution of the static field energy, we can 
calculate the vortex mass (with quantum number n) per 
unit length in the z-direction, i.e.. 



2tt X / dr rS[ 

10 



\Ag 



A" 



(18) 



Note that by definition of the action, the origin of energy 
is zero when the system has no vortex. The vortex mass 
is classified in the following way, 

M„ < 27m (k^ < 1/2 ;Type-I), (19) 
M„ = 27m (k^ = 1/2 ;BPS), (20) 
M„ > 27m (k^ > 1/2 ; Type -II). (21) 

These relations are analytically [13, H^l and numerically 
[sst investigated. FigO and [1] showed the numerically 
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solved K^-dependence of Mn/n, the vortex mass divided 
by the topological number n, and we can see 



< nMi 
> nMi 



< 1/2), 
> 1/2). 



(22) 



These relations indicate that the interaction between vor- 
tices is attractive for <1/2, and repulsive for >1/2. 
These are confirmed by the studies of the vortex-vortex 
interaction analytically in large separation case [36| and 
numerically in arbitrary distance case (ssj . 

The qualitative understanding for the potential prop- 
erty is as follows. The attractive part of the potential 
comes from the reduction of the topological defect of the 
Higgs field condensate. This is because the total energy 
of the system decreases when the domain of the normal 
conducting state decreases. On the other hand, the re- 
pulsive part comes from the magnetic field interaction of 
the vortices. Then, in the case of <1/2 where the co- 
herence length of the Higgs field exceeds the penetration 
depth of the magnetic field, the attractive part becomes 
larger than the repulsive part in the whole region, and the 
vortex- vortex potential becomes attractive. On the other 
hand, in the case of >1/2, the relation is reversed and 
the vortex-vortex potential becomes repulsive. 

In the case of the vortex-anti-vortex case, the poten- 
tial becomes attractive independently from since the 
magnetic field interaction between the vortex and the 
anti-vortex becomes attractive in contrast to the vortex- 
vortex case, and there is no repulsive part in the poten- 
tial. 

Although we consider the isolated vortex under the ide- 
alized situation from the beginning, it is worth mention- 
ing under which situations single vortex appears referring 
to the potential properties explained above. In Type-I 
superconductor, the multi-vortices with total topologi- 
cal number N transfer to the single, giant vortex with 
topological number N due to their attractive interaction. 
After the successive fusion processes, the multi-vortices 
system transfers to the normal conducting state. For the 
finite size superconducting materials under the homoge- 
neous magnetic fields and thermal equilibrium, such fu- 
sion processes lead the difficulty to observe the Type-I 
vortices because the magnetic fluxes successively invade 
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FIG. 3: The plot for dependence of the vortex mass di- 
vided by topological number n, Mn/n for n=l, 2, 3 in the 
rescaled Abelian Higgs model in units of Eq.([2]). 
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FIG. 4: The plot for k dependence of the vortex mass di- 
vided by topological number n, Mn/n for n=l, 2, 3 in the 
rescaled Abelian Higgs model in units of Eq.([2ll. The line 
with circle represents the mass of the n=l vortex calculated 
with Type-II limit (large k'^) approximation in The hor- 
izontal axis represents the log plot of value in the much 
wider range than that in Fig. |3] 



from the interface between the normal and superconduct- 
ing phase. Nevertheless, under the non-equilibrium sit- 
uations or inhomogeneous magnetic fields concentrated 
on the small area, it is expected that the Type-I vor- 
tex or, more generally, the domain of the normal state in 
the superconductor can appear as a transient state in the 
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process such as the dynamical evohition of the thermody- 
namic unstable state into the equilibrium state, i.e., the 
nucleation or spinodal decomposition processes [s^, [11] . 
From the viewpoint of the particle physics, an another 
interesting situation where Type-I vortex appears is the 
case where magnetic monopoles are immersed in super- 
conductors and the isolated vortex emerges between the 
monopole and anti-monopole Q . In this case, no collapse 
with the other vortices occurs and Type-I vortex behaves 
as the stable object due to the topological charge conser- 
vation. 

On the other hand, in Type-II superconductor, the gi- 
ant vortex with topological charge iV fissions to the N 
vortices with charge n=l due to their repulsive nature 
[3^. This repulsive interaction prevents the vortices to 
coalescence and, under the thermal equilibrium, the n—1 
vortices form some stable structure like the Abrikosov 
triangle lattice [l]. Then, in Type-II superconductor, the 
n—1 case is realistic in the thermodynamic equilibrium in 
contrast to Type-I case. However, when the system expe- 
rience a certain quenching [38] such as the rapid change 
of the temperature or density, the n > 2 vortex may be 
generated as the excited states. This is the reason why 
we also consider the Type-II giant vortex in the follow- 
ing. Although they are unstable against the fission into 
small vortices, when is not too large, it may be in- 
teresting to consider the vibrations of the giant vortex 
because such vibrations can occur during the slow devel- 
opment of the fission process due to the small gradient 
of the vortex-vortex potential for the small separation 
of vortices [11]. Moreover, the study of the fluctuation 
modes of the n > 2 vortex may give some insight for the 
early stage of the evolution of the giant vortex into the 
multi- vortices. 



III. AXIAL-SYMMETRIC FLUCTUATIONS 

In this section, using the results of the previous section, 
we analyze axial-symmetric fluctuations in the system 
which has only one vortex in the background. In sub- 
section A, using the static vortex profile, we derive the 
vortex-induced potential for the Higgs and photon field 
fluctuations, and the equation of motion. In subsection 
B, we closely discuss the potential feature to classify the 
type of the excitation modes. 



A. Field equation for the fluctuations 

We investigate the axial-symmetric fluctuations of the 
Higgs field and the related photon fields around the static 
vortex. We expand the fields with small amplitude fluc- 
tuations around the static fleld. 



0(i,z,r) ^0^'(r)+^(i,2,r), 
A0{t,z,r) = A^^ir) +ae{t,z,r). 



Here we only consider the axial-symmetric fluctuations 
and we include t, z-dependence at this stage. It is im- 
portant to notice that, with the Ansatz ^ for the static 
solution, all the other type of fluctuations such as fluctu- 
ations of Aq, Az, Ar and the 0-dependent fluctuations are 
decoupled from the axial-symmetric fluctuations which 
we consider. Hence we can analyze the axial-symmetric 
fluctuations independently from the other type of 
fluctuations and we drop off such fluctuation terms from 
the following Lagrangian. The complete form of the La- 
grangian is given in the Appendix |^ 

Substituting Eq. (P5)l into Eq.® and using the static 
solution, the following expression for the axial-fluctuation 
part of Lagrangian can be derived as 



L 



(2) 
axial 



(3) 
axial 



c 



(4) 

axial ■ 



(24) 



The first term ^l^j^i 



represents ^^Lai = £?^:°i]<^ni,^n 



axial IV^ch^clJ- 

The 1st order terms of fluctuations Lp and ag do not ap- 
pear because we expand around the static solution. The 
2nd, 3rd, and 4th order fluctuation terms take the fol- 
lowing forms; 



1 



1 



C 



(3) _ 

axial 
''-axial- ^ V 



2 

n 



-{drag 



r 



(25) 



Ar 



Lp^ag 



2 2 



20ci¥'a^,(26) 
(27) 



Here V^ip, Vaa, Vipa denote the vortex-induced potential, 
i.e., these interactions are induced by the static vortex 
configuration which can be regarded as the large ampli- 
tude background external field. and Vaa play roles as 
the single particle potential for ip and ag respectively, and 
Vipa is the interaction which mixes the ip and ag. Here- 
after, we will simply call these interactions "potentials" . 
The concrete forms of the potentials are 



V^air) 
Vaair) : 



Ar 



-f 2^2(302^ -1), 



Ar 



(28) 



(23) 



To investigate the fluctuation spectrum, we solve the 
Euler-Lagrange equation for £^2^ neglecting the higher 
order fluctuations. The quantum loop corrections can 
be calculated using the Green's fimctions constructed of 
bases which diagonalize £^2' . However, because of the in- 
homogeneous vortex-background, the Green's functions 
are not functions of the distance between arbitrary two 
points, i.e., G(x,x') ^ G{x — x'), and the calculations of 
the loop corrections are not straightforward. In addition, 
higher order terms of the axial-symmetric fluctuations 
couple with the terms of the non-axial symmetric ones. 
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We reserve the calculation of the loop corrections for the 
future work and, in this work, we safely neglect the quan- 
tum corrections constraining ourselves to the case where 
the perturbative expansion parameters and e^K^ = X 
are sufficiently small. (Recall that the perturbative cor- 
rections are calculated with expansion of e*'^ = e**^ 
See also Eq.Q.) These constraints are not so strict be- 
cause the Higgs-photon coupling e is usually small and 
we mainly consider ^ 1/2 cases, where e^K^ = A is 
also sufficiently small. 

For later convenience to discuss the effect of the po- 
tential, we change the variables 



These changes of variables enable us to considerably sim- 
plify the equation of motion as will be seen. The Euler- 
Lagrange equation for 2nd order Lagrangian is obtained 
through the variation. 



f dtdzdr2TrrC'^^^ = — [ dtdzdr2TirC^'^^ = 0. (30) 
oae J 



5(p 



<f{t,r,z) 
ae{t,r,z) 



- .-1/2 



5/V2 



(29) Then we obtain the following coupled equation 



if{t,z,r) 
ag{t,z,r) 



-dl 



2v^0d(n-A^i) _q2 I _3_ i o^2 



3 
4^2 



20^ 



'f{t,z,r) 
ag{t,z,r) 



-a2 + y„„(r) V,,(}{r) 



(fit,z,r) 
ag{t,z,r) 



(31) 
(32) 



Here Va/sir) — Vpa{r) (This is the consequence of the 
change of the variables in Eq. ([29|) ). Note that the t, z 
dependence of (p and dg is trivial because the background 
static vortex configuration is t, z-independent and t, z de- 
pendence does not appear in the potential. Then the par- 
tial differential equation can be decoupled into the total 
differential equations for t, z, r, respectively. 

The system has the translation invariancc in both z- 
direction and time direction, and then the solution takes 
the following plane wave solution, 

( f{t,r,z) \^_,^^,_,^^)/ «(r)\ 
\^ ag{t,r,z) J \ I3{r) J 

Then the total differential equation in the radial direction 
for a{r) and /3(r) is obtained as 



/3(r) is expressed as 

/ dr{a,,p,){ak,f3kf = S,k / dr [a,{rf + (3,(rf] 
Jo Jo 

= Sjk- (35) 

Here we give a brief explanation for the qualitative pic- 
ture for the excitation modes. The shape of the radial 
direction is given by aj{r) and I3j{r). This wave function 
propagates in the z-direction with the frequency 0Jj,k^ 
and the wave vector fc^. Therefore the root of the eigen- 
value m?j of the radial wave functions aj [r] and 13 j (r) play 
roles of the mass for the propagation in the z-direction. 
Then we call this type of new mode as the peristaltic 
mode of the vortex as illustrated in FiglTJ 



(34) 



kl- Note that the Hermiticity of the 



where m| 

matrix in the RHS in Eq. p4|) ensures the real value of 
and our later analysis shows m| > 0, i.e., the stability of 
the static profile against the fluctuations. The root of the 



eigenvalue for the Eq. ((34| 



plays a role of the mass for 



the propagating mode in the z-direction. The Hermiticity 
of Eq. (|M|) also enables us to take the functions a(r) and 
(3{r) real, and the orthonormal condition for the a(r) and 



B. The feature of the potential 

We show in FigO the single particle potentials Vaa , 
Vaf3 and Vf3f3 around the vortex for various values of k. 
We can deduce the considerable information of the exci- 
tation modes from the examination of these potentials. 
To get the qualitative picture about the wave functions 
of fluctuations, it is convenient to divide the region of the 
coordinate space r into three regions, central (0 ^ 0.5), 
intermediate (0.5 4), and asymptotic (4 ^) region in 
the rescaled Abelian Higgs model. 

In the central region, when r — > 0, the asymptotic be- 
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FIG. 5: The plot of the vortex-tnduced potential Vaa{r),Vaj3(r),V()i3{r) as the single particle potential for fluctuations a{r) 
and /3(r) around single vortex, defined in Ea. (l33p . The Ginzburg-Landau parameter takes k^=0.1, 0.5, 1.0, and the topological 
number takes n=l. 



havior of the potential is 

- 1 3 

Ka;/3 = Vfja — * (negative const.) x r"~"^ . (36) 

From these asymptotic behavior, it is expected that the 
fluctuations a{r) and (3{r) around the vortex center are 
suppressed by the strong repulsive interaction propor- 
tional 1 /r^ and their distributions are shifted to the out- 
side. 

On the other hand, in the intermediate region, the 
strong repulsive part of Vaa and Vpp becomes weak and 
the attractive pocket emerges. The mixing term Vap also 
plays an important role to enhance the mixing between 
the Higgs and photon field fluctuations, or a and (3. Since 
the Vap always takes the negative value, the excitation 
energy is reduced when the relative sign of the a and (3 
fluctuations is positive. To avoid the confusion, we stress 
that the minus sign of Vap does not mean attractive inter- 
action. Even if Vap takes the positive value, the energy is 
reduced when the relative sign of a and (3 fluctuations is 
negative. This is also recognized from the fact that the 
replacement (a,/3) {a,—f3) changes only the sign of 
the off-diagonal part of Ea. ([M)) . In treatment up to the 
second order of the fluctuations, only the absolute value 
of Vap is relevant to the mixing of a and /?. All these 
features of the potentials Vaa,Vpp, Vap make low energy 
modes of the fluctuations localized in the intermediate 
region. 

Finally we discuss the asymptotic region. When r 
oo, these potentials asymptotically behave as 

ya„^4K2, Vpp^2, Vo,p = Vpo,-^0. (37) 

Note that the asymptotic behaviors do not depend on 
the topological charge n. Since Vap goes to zero, the 
correlation between a and (3 vanishes and these fluctua- 
tions behave as independent modes. Then we have only 
to consider the one-particle problem in this region, i.e.. 



we treat the asymptotic forms of Eq. ([M]) as 

(38) 

(-^ + 2)/3,(r)^m|/3,(r). (39) 

The values and 2 in the LHS play roles of the (square 
of) continuum thresholds of the a and (3 respectively. 
These thresholds are equivalent to the masses of the fluc- 
tuations in the system with no vortex, i.e., 2k corre- 
sponds to the mass of the Higgs field and V2 is the photon 
mass generated through the Anderson-Higgs mechanism 
[3S . J:0] . Then, the threshold energy for the continuum 
states is expressed as 

TOth = min(2K, \/2). (40) 

When m?j is smaller (larger) than the thresholds, the fluc- 
tuations behave like bound (continuum) states. 

It is also worth mentioning that Vaa and Vpp take the 
same value in the asymptotic region and the continuum 
thresholds are the same both for a and (3 for the case 
of K^=l/2, where BPS saturation occurs. This note is 
sometimes useful to roughly classify the type of the fluc- 
tuation spectrum in the low-energy region in terms of 
the type of superconductor, as is explained below. It is 
expected that, for Type-I case (k^ <1/2), the threshold 
for the Higgs field is lower than that for the photon field, 
and the Higgs field fluctuation dominates over the low 
energy physics. For Type-II case {k^ >1/2), the rela- 
tion is reversed, the photon field fluctuation dominates. 
We also expect that the collective properties induced by 
the mixing of the Higgs and photon field fluctuations are 
most pronounced when the BPS saturation occurs. This 
is because the strength of both fluctuations are the same 
order and they affect each other. In the next section, 
we will see the above mechanism is crucial to emerge the 
characteristic discrete pole around = 1/2. 



9 



1 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



;n=l) ot-thres. 

P-thres. -o- 



j \ 1 r 

"fcWWVV VV V V V^' 
ililll i I I i ▲ A AiHU 



AMfWAAA A/ 



□ jn imiii II ii mrmmrnf ] 

><»»0<XX X X X ><mXX8»iX>0<XXX< XX XXX XK ■> 

^H-++++ + + ++ l imn I II I 1 1 I +H-H--I+ +++ + 



_i_ 



_i_ 



0.5 1 1.5 



2.5 3 3.5 



m 







n 1 1 — 

7VV VV VW ^WMW 



FIG. 6: K, -dependence of rrij for n—1 case and the square 
of thresholds for a and /3 fields, and 2, respectively. 
Around k^=1/2 (BPS value), as a collective mode of the Higgs 
and photon field fluctuations, there appear discrete poles in 
lower-energy region below (square of) the continuum thresh- 
old mtii = min(4K^, 2). 
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FIG. 7: The same plot as Fig|6] except the replacement of 
rrij with mf, i.e., calculated in the case that the channel 
coupling between Higgs and photon field fluctuations is ar- 
tificially dropped off. Without channel coupling, the bound- 
state with ~ 1/2 becomes the shallow bound state due to 
the lack of the interplay between the Higgs and photon field 
fluctuations. 



IV. APPEARANCE OF THE DISCRETE POLE 
AS THE PERISTALTIC MODE OF THE VORTEX 



In this section, we summarize the numerical results in 
both Type-I and Type-II superconductors. The follow- 
ing results are obtained by solving the coupled channel 
equations ([M]) numerically with the lattice spacing 0.02 
and the lattice size 40.0 in our rescaled unit. We impose 
the boundary condition on the wave function a{r) and 
/9(r) to be zero value at the box boundary. In the fol- 
lowing, we consider kz = case. The energy and radial 
wave functions shown in this section are translated into 
those in the general kz ^ case with the replacements 



E = rrij ^ {m'j + kiy/"^ and multiplying e^^'^ to the 
radial wave functions, respectively. 



A. Analysis for the excitation spectrum 

In this subsection, we show the numerical results for 
n=l case in both Type-I and Type-II superconductors. 
In FiglSl the K^-dependence of m| for n=l case and the 
thresholds for a and (3 fields are plotted. Shown in Fig[7] 
are the same quantities as Fig|Sl but without channel 
coupling, which is artificially dropped off. 

With smaller k^(< 1/2), the threshold for a is shifted 
lower than threshold for /?, and a can excite with smaller 
energy and the amplitude of a is enhanced compared 
to that of (3 field. Then the lower excitation modes are 
mainly characterized by the a field fluctuation. Under 
the normalization condition ([55]) . there appears a large 
excess of the amplitude of a field compared to that of 
f3 field, and this mismatch between a and (3 leads the 
suppression of the mixing effect. Then the properties of 
the lower excitations asymptotically approach to those 
calculated without channel coupling. The property of 
the discrete pole can be interpreted as the bound state 
of the static vortex and the fluctuating Higgs field. 

On the other hand, with larger k^(> 1/2), the enhance- 
ment of the threshold energy for a leads the suppression 
of the a field excitation, and simultaneously enhances 
the (3 field ratio in the lower excitation modes. Then the 
increasing of value leads the reduction of the mixing 
correlation energy which is indispensable for the appear- 
ance of the discrete pole. Eventually, at ~ 0.8, the 
discrete pole is buried in the continuum state spectrum of 
the fluctuating photon field and no discrete pole appears 
for > 0.8. This indicates that such a pole is difficult 
to be identified in most Type-II superconductors. 

Finally we consider the most interesting case, ~ 1/2 
case, i.e., near BPS saturation. In this case, the mixing 
effect becomes larger, and, as a result, the correlation 
between a and (3 considerably lowers the excitation en- 
ergy. This is clearly seen from the comparison of Fig [6] 
with FiglTl In Table |T1 we summarize the results for 
the lowest excitation energy m^^o f^nd their relations to 
the threshold energy mth = min{2K, v^}, the lowest ex- 
citation energy without channel coupling, TOq, and the 
energy differences among them. The quantity mth ~ "^o 
represents the binding energy and specifies the spatial 
extension of the wave functions. The quantity TOq — toq 
is, roughly speaking, the correlation energy induced by 
the mixing effect between Higgs and photon field fluctu- 
ations. As expected, the maximum mixing correlation is 
realized in the BPS case (k^=1/2). 
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FIG. 8: The behavior of the low-lying mode of Higgs field fluctuation ip{r) = r^^'^a{r), the photon fluctuation ae(r) = 
y/2r^^'^ P{r) and the magnetic field fluctuation fez(r) = ^-^{rag) for n=l vortex. For — 0.1, 0.5 cases, we plot the lowest 
excitation mode with the radial mass mo smaller than the threshold mth. For = 1.0 case, where the discrete pole does not 
appear, we plot the typical low-lying excitation mode with the radial mass rrij, satisfying 2 < < 4k^. 



TABLE I: dependence of the various energy in the rescaled 
unit. The physical value can be obtained by multiplying v/e. 
Tith, mo, m'o represent the threshold energy, the lowest ex- 
citation energy and the lowest excitation energy without the 
channel coupling, respectively. The negative sign of mo — mth 
indicates the bound state. The quantity mo — mo roughly 
characterizes the correlation energy of the Higgs-photon field 
mixing. 





mth 


mo 


mo — mth 


m'o 


mo — mo 


0.3 


1.095 


1.003 


-0.092 


1.068 


-0.065 


0.4 


1.265 


1.138 


-0.127 


1.229 


-0.091 


0.5 


1.414 


1.247 


-0.167 


1.371 


-0.124 


0.6 


1.414 


1.333 


-0.081 


1.414 


-0.081 


0.7 


1.414 


1.394 


-0.020 


1.414 


-0.020 
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m 
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B. The spatial behavior of the fluctuations and the 
character change of the excitation modes with k 



FIG. 9: The /t^-dependence of = dr a(r)^ 

dr Pir)"^ and A - B = dr a{r)f3{r) for the lowest ei 
tation mode in the case of n = 1. 



The features explained above can be explicitly seen 
from the behavior of the radial wave function of the Higgs 
and the photon fields, i.e.. 



(^(r) = r^/'^a{r), ae{r) = V2r^/'^ f3{r). 



(41) 



For the argument of the radial part, we hereafter drop 
off the trivial factor e^*'^*+'^^^ appearing in the relations 
(|29)) and p3|). In fact, ip{r) and ae{r) denote the radial 
part of the Higgs and the photon fluctuations in Eq. (P5|) . 

In FigjS] we show the wave functions (p{r) and ae{r) 
of the low-lying modes and the corresponding magnetic 



fluctuations 



hzir) = --^{raeir)), 
r dr 



(42) 



for k'^= 0.1, 0.5, 1.0. For =o.1 and 0.5, we plot the 
lowest eigenfunction with radial mass smaller than the 
threshold mth- On the other hand, for — l.O, there is 
no discrete pole because the radial mass mj is larger than 
the threshold mth = V^, i.e., > 2. Then we plot the 
typical wave function of the low lying mode whose radial 
mass satisfying 2 < nij < Ak^. 
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The behavior of the fluctuation fields coincides to our 
previous expectations. For k^—0.1, the Higgs field fiuc- 
tuation considerably exceeds to the photon fluctuation. 
Although both of fluctuation fields have the smaller exci- 
tation energy than the threshold mth and this leads local- 
ization of them around the vortex core, the distribution 
of the Higgs field fiuctuation spreads broader than that 
of the photon field fluctuation. This is because the low- 
est excitation energy toq is close to the threshold for the 
Higgs field fluctuations, 2k, and the lowest mode is char- 
acterized as the shallow bound-state between the static 
vortex and the Higgs field fiuctuation. 

When we change to the k^=0.5, the mixing correla- 
tion induces the corporative behavior of the Higgs and 
photon field fluctuations and enlarges the binding energy 
mth — 'niQ. As a result, the lowest excitation mode be- 
comes the deeper bound-state, and then, the Higgs field 
localizes closer to the core, which affects the magnetic 
field vibration more strongly. 

For K^=1.0, because of 2 < rnj < Ak^ for low-lying 
modes, only Higgs field fluctuation is localized and the 
photon field ao{r) oscillates asymptotically (r — > oo) as 



ae{r) 



1/2 



(43) 



leading to the asymptotic behavior of the magnetic field 
hz(r) like 



r dr 



(44) 



with = m'j — 2. (More precisely, ag{r) and hz{r) are 
expressed as the linear combinations of the cos(fcrr)/r^/^ 
and sin(fcj,r)/r^/^ because (p{r) and ag(r) are real func- 
tions.) 



Here we comment on the total magnetic fiux. In the 
static case, the total magnetic fiux takes the quantized 
value 27rn, as seen in Eq.([T7]). On the other hand, the 
photon fiuctuation generally induces a time-dependent 
magnetic flux, whose proper time-average vanishes due 
to the factor e~*"*. The contribution from the radial 
part of the magnetic field fluctuation is 



(5$= / dr 2'!rrhz{r) = / dr 2Tr-^{rae{r)) 



= 27r[ra0(r)]^ 



(45) 



For r — > 0, ae(r) asymptotically approaches to zero. On 
the other hand, for r — > oo, 



rae{r) 



(discrete pole), ^-^g-j 
(continuum states). 



where C, D represent some constants. Then, for the 
discrete pole, (5$ becomes zero and thus the total mag- 
netic fiux is conserved in the arbitrarily cross section 
perpendicular to the vortex line. In contrast, for the 
continuum states, the total magnetic flux becomes sen- 
sitive to the situation far from the vortex core because 
of the long tail behavior of hz{r). For the continuum 
state, the radial photon fluctuation ae{r) asymptotically 
behaves as a two-dimensional spherical wave in the ho- 
mogeneous no- vortex vacuum, and hence its asymptotic 
part can be regarded as a non-vortex-origin fluctuation. 
In other words, the continuum states can be regarded as 
the vortex plus a "scattering wave" . 

To discuss the character change of the excitation modes 
in wider range of k^, we deflne the functions 



dr a(r)^. 



dr I3{rf 



(47) 
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Notice that dr [a{r)^ + /?(r)^] = 1 in the normahza- 
tion condition. fSee Eg. ([551).) We also define the function 
to estimate the K^-dependence of the collective properties 
and their mixing ratio for the lowest eigenmode. One of 
the good indicators to estimate the degree of the collec- 
tive nature is the function 



A-B = 



dr a{r)(3{r). 



(48) 



When the oscillation is sympathetic, this integral takes 
the large value. In Figl9]we plot the behavior of the func- 
tions A^, and A ■ B for various . As expected, the 
ration of a(r), or Higgs field decreases compared to that 
of (i{r), or photon fluctuation with increasing because 
the threshold for the Higgs field fluctuation increases. 
Eventually, around ~ 0.8, the discrete pole buries in 
the continuum, and the lowest mode becomes the con- 
tinuum state of the photon field. At the same time, the 
correlation between the Higgs and photon field fiuctua- 
tions becomes very small, as seen in the lower panel in 
FigH 



C. The peristaltic modes for n > 2 cases 

In this subsection, we show the results for the peri- 
staltic modes of the n > 2 vortex, which is realistic in 
the case of < 1/2. We emphasize the difference of 
the static vortex between n=l vortex and n=2, 3 vortex, 
and then discuss their effects on the potential behavior 
and the resulting excitation spectrum. In this subsec- 
tion, when we argue n-dependence of the classical profile 
and the potential, value is fixed to BPS value, i.e., 
k2 = 1/2. 

Shown in Fig[Tn]are the n=l, 2, 3 static vortex profiles 
and the corresponding vortex-induced potentials Vaa, 
Vpp, Vap for the fluctuations a(r) and /3(r), respectively. 
Here we qualitatively explain the changes of the static 
vortex profiles in terms of n. Since the total penetrating 
fiux of the magnetic field is proportional to the topologi- 
cal number n, the magnetic field around the vortex core 
increases with increasing n. Then the enlarged magnetic 
flux pushes the Higgs field outward or reduce the ex- 
pectation value of the Higgs fleld. Roughly, this process 
spreads the normal phase region near the vortex core and 
leads the outward shift of the interface between normal 
and superconducting phase. According to the outward 
shift of the surface region, where the kinetic energy of 
the Higgs field is enhanced, the maximum of the local 
energy density also shifts outward, as seen in the upper 
panel of FigfTOl 

Next we present the qualitative explanation for the sin- 
gle particle potential for the fluctuation fields. First we 
give the explanation for the change of Vaa for the Higgs 
field fluctuation. Since the area of the magnetic flux is 
spread with increasing n, the intersection surface of the 
static magnetic and the Higgs fields moves outward. This 
leads the outward shift of the potential Vaa- 



The changes of V3/3 can be interpreted in somewhat in- 
direct way as follows. The fact read off from Fig.fTUlis the 
depth enhancement of the attractive pocket in Vpp with 
increasing n. This enables the photon field fluctuations 
to enhance and localize around the surface region more 
easily for larger n cases. The localization of the photon 
field fluctuation means somewhat rapid spatial change of 
the wave function ag{r), and leads a large fluctuation of 
the magnetic field hz{r) = ^^(rae). Then we can say 
that the depth enhancement of the attractive pocket in 
V/3/3 for large n leads the large fiuctuation of the mag- 
netic field. The large magnetic fluctuation reflects the 
reduction of the Higgs field condensate or the weakened 
Meissner effect. In fact, the depth enhancement of the 
attractive pocket is closely related to the outward shift 
of the static Higgs field distribution. 

Finally we comment on the potential Vap, which in- 
duces the mixing correlation between the Higgs and the 
photon field fluctuations. As seen in Fig. [TOl with in- 
creasing n, the maximum of {Vapl shifts outward. The 
physical meaning of the change in this potential is natu- 
rally explained as the outward shift of the interface be- 
tween the static Higgs field and magnetic field because, 
at the interface, the fiuctuations of the Higgs and photon 
field can frequently interplay. 

With increasing n, the resulting changes in the static 
configurations and the vortex-induced potential lead the 
change of the excitation energy of the peristaltic modes. 
Now we see the excitation spectra of the giant vortex 
(n — 2,3) and their dependence on the value. The 
numerical results for n = 2, 3 cases are shown in Figs. lTTl 
1121 respectively. For larger n, the excitation energy for 
the lowest mode becomes much smaller as expected from 
the arguments in the previous paragraphs. As a conse- 
quence of the depth enhancement of the attractive pocket 
in the discrete pole appears up to considerably large 
value in the Type-II region. 

In such large k^(> 1/2) cases, of course, the fission pro- 
cesses of the giant vortex to many n=l vortices become 
important. According to the analysis for the instability 
of the static vortex [23, |4l| and the interaction between 
two vortices 33| , there exists no potential barrier for the 
fission from single Type-II giant vortex to many n=l vor- 
tices and the giant vortex can fission through the non- 
axial symmetric modes. However, in spite of this observa- 
tion, we have a possibility that axial-symmetric discrete 
poles play important roles if the lifetime (decay width) 
of the giant vortex is very large compared to the typi- 
cal time scale (excitation energy) of the axial-symmetric 
fiuctuations. The decay rate can be calculated through 
the analysis for the non-axial-symmetric modes and the 
fiuctuations of AQ,Ar, which are not considered in this 
work. We expect that, for not so far beyond the BPS 
value, — 1/2, the lifetime is large because the force to 
break up the giant vortex is small. Then, if the lifetime is 
confirmed to be large enough, we can immediately apply 
our results without further consideration about the corre- 
lation between the axial-symmetric and the other modes 
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because these modes behave independently as far as the 
distortion is sufficiently small, as noticed in Sec.IIIA. 

On the other hand, for < 1/2 cases, since the gi- 
ant vortex is stable and the axial-symmetric modes are 
independent from the other modes not considered in this 
work, the strong bound-state-like modes more easily ap- 
pear than the case for > 1/2 because of non-existence 
of the fission processes. In addition, according to the 
negative surface tension, vortex with k? < 1/2 tends 
to reduce the surface area and then the axial-symmetric 
fluctuation becomes increasingly prominent with smaller 
K^. Therefore, we expect that the discrete pole shown 
in this section shares the considerable relevance to the 
low-energy behavior of the giant vortex with < 1/2. 



V. SUMMARY AND CONCLUDING REMARKS 

In this paper, as a new-type mode beyond the stringy 
one, we have studied the axial-symmetric peristaltic 
modes of the vortex, peculiar to the extension of the 
flux-tube. We have investigated the peristaltic modes 
both in Type-I and Type-II superconductors for n —I, 2, 
3 and = 0.1 ~ 1.0. We have treated single vortex in 
the situation where Type-I vortex can be realized with- 
out collapses with the other vortices and behaves as the 
stable object. 

Here we summarize our main results about the fluc- 
tuation modes. With < 1/2, the lowest excitation 
modes are characterized by the Higgs field fluctuation 
binding with the static vortex. On the other hand, with 
> 1/2, the mass of the Higgs field fluctuation becomes 
heavier and the photon fluctuation modes are relevant. 
For = 1/2 (BPS cases), there appears the characteris- 
tic discrete pole. In this case, the Higgs and the photon 
field fiuctuations have the same masses, and then both 
fiuctuations play important roles and their corporative 
behavior makes the lowest excitation softer one. This 
mechanism is the main origin of the characteristic dis- 
crete pole. Since the energy separation between this dis- 
crete pole and the continuum threshold is large around 
~ 1/2, the peristaltic mode may appear most clearly 
in the superconductor with ~ 1/2. 

We have also investigated the fluctuation modes of the 
giant vortex with the topological number n > 2. It is 
found that, with larger n, the photon field fluctuation 
becomes more important for the lowest eigenmodc be- 
cause the distribution of the static Higgs field is shifted 
outward and the photon fluctuation can excite easily 
around the core. We expect that, in < 1/2 cases, the 
axial-symmetric excitation mode with the bound-state- 
like pole plays an important role in the low-energy dy- 
namics of the vortex. On the other hand, for > 1/2 
cases, although the equations for axial-symmetric ex- 
citation modes are decoupled from those of the other 
modes, it is necessarily not only to investigate axial- 
symmetric modes but also to consider whether these ex- 
citation modes can be realized or not during the lifetime 
of the giant vortex. To estimate the lifetime, we have 
to treat the non-axial-symmetric modes and the fluctua- 
tions of Ao,Ar, which drive the fission of the single giant 
vortex to many n=l vortices. We expect that, for 
value not so far from the BPS value, k^=1/2, the life 
time is long enough and we can explore the possibility of 
the appearance of the discrete pole. 

Finally, we end with a future perspective toward the 
further investigation for the dynamics of the vortices (or 
color-electric flux tubes) regarding them as dynamical 
degrees of freedom. 

Although we restrict ourselves to the axial-symmetric 
static vortex profile and fluctuations of 4> and Ag, we 
will investigate the more general cases without using the 
axial-symmetric Ansatz for the static profile and clas- 
sify the large class of the fluctuation modes, i.e., {(j), Aq, 
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Ag, Ar) in X = gauge, or equivalently, {(j), x, Ag, Ar) 
in Aq = gauge, which are not fully considered in this 
work. For this purposes, we will calculate the static vor- 
tex profile and excitation spectrum of the fluctuations in 
two-dimensional analysis in near future. 

For Type-I case, we expect that, in low energy, the 
relevant excitation mode around single vortex is axial- 
symmetric since the Type-I vortex tends to reduce the 
surface area owing to the positive surface tension. Also 
from the point of view on the kinetic energy, the axial- 
symmetric excitation is favorable. On the other hand, the 
arguments for the Type-II case are more subtle. Because 
Type-II vortex tends to increase the surface area accord- 
ing to the negative surface tension, it may be important 
to compare the surface energy to the kinetic energy of 
the distortion. 

For n > 2 giant vortices, interesting subjects are re- 
mained. One of them is the calculation of the lifetime of 
single giant vortex for > 1/2 and the comparison to 
the excitation modes considered in this work. The life- 
time of the giant vortex is closely related to the dynamics 
of the multi- vortices in non-equilibrium systems. As well 
as the information for the vortex-vortex potential, the 
information for the lifetime of single giant vortex may 
provide the good building block to understand the inter- 
esting but complicated phenomena such as the transition 
of the system from the thermodynamically unstable state 
to the stable one. In connection with them, we are also 
interested in the fluctuations between multi-vortices and 
their evolution through the vortex- vortex scattering [i^ . 
As a successive study, we would like to investigate the 
results of the full two-dimensional analysis for excitation 
modes of the vortex in both single vortex and the multi- 
vortices system 

In application to QCD, the dynamics of the color- 
electric flux-tubes, or gluonic excitations as nonpertur- 
bative modes are also important to discuss the hadrons, 
which are composed of quarks and color vortex between 
quarks. In particular, the analysis for the multi-color- 
flux-tube system can give an good insight to understand 
the color confinement- deconfinement transition in QCD, 
and properties of the hot and dense matter of hadrons 
or quarks. It is also interesting to consider heavy ion 
collisions and the successive occuring expansion of the 
hot quark matter in close analogy with the evolution of 
the early Universe. The distribution of the topological 
defects may affect collective properties of quark gluon 
plasma and the resulting particle productions. 

In this work, we have studied the dynamics of single 



vortex as one example of the dynamics of the topological 
objects. The treatments of topological objects as essen- 
tial degrees of freedom are quite general and important 
to understand the crucial aspects of the nonlinear field 
theory, which sometimes can not be reached with the per- 
turbative treatments for the nonlinear terms. We would 
like to extend our analysis to the other topological objects 
such as color flux-tubes, instantons, magnetic monopoles 
in gauge theories, or branes in the string theory and 
their many-body dynamics. 
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Note added: After this work is completed, we notice a 
similar study using a different gauge where ghosts appear 
p5| . Their numerical result on the physical modes is 
consistent with ours. 



APPENDIX A: THE DECOUPLING OF THE 
AXIAL-SYMMETRIC FLUCTUATION MODES 

We give the complete form of the action in the cylindri- 
cal coordinates (r, 0, z). We use the polar-decomposition 
of Eq. dn]) for the Higgs field 

i^{x) = (t){t, r, e, 2)e™«+^x(*^'--»^^) (-TT < X < tt), (A1) 

and fix the gauge to x = 0. Then we obtain the action 
as 

S^^"^ = jdtdrdOdz r£X=°, (A2) 

where 



1 



1 



{n-rAg)' 



4' + {Al-Al 



(A3) 
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Here the electric and magnetic fields are described as 



E,=dtAr-drAo, Eg ^ dtAg - -dgAo, ^ dtA, ~ O^Aq, 

r 

Br^-deA,-d,Ag, Bg = d,Ar-drA,, = -deirAr) - -dgAr 



We decompose the field into the static part and fluctuation part, 

A,{t,r,e,z)^ Af {r,e,z)^aj{t,r,e,z), 
where j = (0, r, 9, z). For the static part, we impose the Ansatz ([8]), or equivalently, 

q^-\r, e, z) = 0^'(r), Af{r, 0, z) - Af{r), < = Af = Af = 0. 
Putting these forms into the Lagrangian (|A3|) . we obtain 

Ag] = C^' + £(2) + £(3) + £(4)^ 



(A4) 
(A5) 

(A6) 
(A7) 
(A8) 

(A9) 

(AlO) 

(All) 
(A12) 

The Lagrangian for the static part, i.e., /C^' is related to static energy pT|) as £ = — £'^'. The fluctuations of the electric 
and magnetic fields, i.e., (6^,6^,62) and {bj.,bg,bz) are expressed as the Ea. (IA4l) with the substitution Aj Uj. The 
linearized Euler-Lagrange equations for the fluctuations ip, ag, ao, ar, Qz are obtained as 



where 



2 r 



{n - rAf) 



cl\2 



ipag, 



£(3) ^ -4k20<=1(^3 ^ 



2^c. 3 , 2(n-rO 



/-(4) 2,' 2 2 2 2\ 2 4 



Q^2 r dr dr' 96*2 
r _ \d_, d_. 
dt"^ r dr^ dr' 



021 

dz"^ 

r dr^ dr' r'^ dO"^ dz"^ 
dt^ ^2 ^5/2 ^^2 



ag = — 



in - rAf) 



\\2 



+ 2k2(30'^i _ 1)2 
H''\n- rAf f 



2(f>''\n~rAf) 
V H ae, 



1 5 



-Lp- 



ao = - 



20cl«r 



9 r 1 9ae 1 d{raj.) da^ 

2 9ae _a 
r2 96' 9r 



9ao 1 d{rar) da 
~dt 



dr 



dz 



52 



19,9 



9t2 r dr^ dr 



J_d^ 
r2 96*2 



-20^10^ 



_9 

dz 



dz 

dao 1 dag daz 
dt r do dz 
dao 1 dag 1 d(rar) 
~dt ^ r~d9 ^ 

I 



dr 



(A13) 
, (A14) 
(A15) 
(A16) 
(A17) 



It is important to notice that we can decouple the equa- 
tions for the axial-symmetric fluctuations of (p and ag 
expanding the fluctuation modes with the Legendre poly- 
nomial. If we restrict ourselves to the / = modes analy- 
sis, we can neglect angular-derivative terms in Eq. (jA13p - 
(|A17[) . Therefore, we can treat the fluctuations (p and 



ag independently from fluctuations aQ,ar,az, because ip 
and ag are related to oq, a,-, only through the angular- 
derivative terms. As a result, to investigate the axial- 
symmetric fluctuations of ip and ag , we have only to treat 
the Lagrangian (P5)) as mentioned in Scc.IIIA. 



[1] A. A. Abrikosov, Soviet Phys. JETP 5, 1174 (1957). 



[2] P. A. M. Dirac, Proc. Roy. Soc. A133, 60 (1934), 



16 



[3] G. 't Hooft, Nucl. Phys. B79, 276 (1974), A. M. 

Polyakov, JETP Lett 20, 194 (1974). 
[4] A. A. Bolavin, A. M. Polyakov, A. S. Schwartz, and Y. 

S. Tyupkin, Phys. Lett. 59B, 85 (1975). 
[5] R. Rajaraman, Solitons and Instantons, (North-Holland, 

1987). 

[6] N. Manton and P. Sutcliffe, Topological Solitons, (Cam- 
bridge University Press, Cambridge, 2004) and references 
therein. 

[7] H. B. Nielsen and P. Olcscn, Nucl. Phys. B61, 45 (1973). 
[8] Y. Nambu, Phys. Rev. DIG, 4262 (1974). 
[9] G. 't Hooft, Nucl. Phys. B190, 455 (1981). 

[10] S. Mandelstam, Phys. Rept. C23, 245 (1976). 

[11] Color Confinement and Hadrons in Quantum Chromo- 
dynamics edited by H. Suganuma, N. Ishii, M. Oka, H. 
Enyo, T. Hatsuda, T. Kunihiro, and K. Yazaki (World 
Scientific, Singapore, 2004) 

[12] H. J. Rothe, Lattice Gauge Theories: An Introduction, 
3rd edition (World Scientific, Singapore, 2005) 

[13] M. Creutz, Phys. Rev. Lett. 43, 553 (1979); 43, 890 
(1979); Phys. Rev. D 21, 2308 (1980). 

[14] T.T. Takahashi, H. Matsufuru, Y. Ncmoto, and H. Sug- 
anuma, Phys. Rev. Lett. 86, 18 (2001); T.T. Takahashi, 
H. Suganuma, Y. Nemoto, and H. Matsufuru, Phys. Rev. 
D 65, 114509 (2002). 

[15] H. Ichie, V. Bornyakov, T. Streuer, and G. Schierholz, 
Nucl. Phys. A721, C899 (2003). 

[16] A. Vilcnkin and E. P. S. Shellard, Cosmic Strings and 
Other Topological Defects, (Cambiridge University Press, 
Cambridge, 1994) and references therein. 

[17] T. W. B. Kibble, J. Phys. A 9, 1387 (1976). 

[18] W. H. Zurek, Nature (London) 317, 505 (1985); Phys. 
Rep. 276, 177 (1996). 

[19] See, for example, M. Tinkham, Introduction to Supercon- 
ductivity (Krieger, Malabar, 1985). 

[20] M. Hindmarsh and A. Rajantie, Phys. Rev. Lett. 85, 
4660 (2000). 

[21] G. J. Stephens, L. M. A. Bettencourt, and W. H. Zurek, 

Phys. Rev. Lett. 88, 137004 (2002). 
[22] For instance, J. Polchinski, String Theory, (Cambridge 

University Press, Cambridge, 1998). 
[23] K. Nawa, H. Suganuma and T. Kojo, Phys. Rev. D85, 



086003 (2007), and its references. 
[24] M. Luschcr, G. Munstcr and P. Wcisz, Nucl. Phys. B180, 
1 (1981). 

[25] P. Orland, Nucl. Phys. B428, 221 (1994) and its refer- 
ences. 

[26] E.T. Akhmedov, M.N. Chernodub, M.I. Polikarpov and 
M.A. Zubkov, Phys. Rev. D53, 2087 (1996) and its ref- 
erences. 

[27] E. B. Bogomol'nyi, Sov. J. Nucl. Phys. 24, 449 (1976). 
[28] M. K. Prasad and C. M. Sommerfield, Phys. Rev. Lett. 
35, 760 (1975). 

[29] H. Suganuma, S. Sasaki, and H. Toki, Nucl. Phys. B435, 
207 (1995). 

[30] H. Suganuma, S. Sasaki, H. Toki, and H. Ichie, Prog. 

Theor. Phys. Suppl. 120, 57 (1995). 
[31] K. J. Juge, J. Kuti, and C. J. Morningstar, Nucl. Phys. 

(Proc. Suppl.) B63, 326 (1998); Phys. Rev. Lett. 90, 

161601 (2003). 

[32] T.T. Takahashi and H. Suganuma, Phys. Rev. Lett. 90, 
182001 (2003); Phys. Rev. D 70, 074506 (2004). 

[33] L. Jacobs and C. Rebbi, Phys. Rev. B19, 4486 (1979). 

[34] H. J. de Vega and F. A. Schaposnik, Phys. Rev. D14, 
1100 (1976). 

[35] P. Tholfsen and H. Meissner, Phys. Rev. 169, 413 (1968). 

[36] L. Kramer, Phys. Rev. B 3, 3821 (1971). 

[37] H. Frahm, S. Ullah and A. T. Dorsey, Phys. Rev. Lett. 

66, 3067 (1991). 
[38] F. Liu, M. Mondello and N. Goldcnfcld, Phys. Rev. Lett. 

66, 3071 (1991). 
[39] P. W. Anderson, Phys. Rev. 130, 439 (1963). 
[40] P. W. Higgs, Phys. Rev. Lett. 13, 508 (1964). 
[41] S. Gustafson and I. M. Sigal, Commun. Math. Phys. 212, 

257 (2000). 

[42] E. Myers, C. R. Rebbi, and R. Strilka, Phys. Rev. D45, 
1355 (1992) 

[43] T. Kojo, H. Suganuma, and K. Tsumura, in preparation. 
[44] For recent review, M. Eto, Y. Isozumi, M. Nitta, K. 

Ohashi, and N. Sakai, J. Phys. A39, R315 (2006). 
[45] M. Goodband and M. Hindmarsh, Phys. Rev. D52, 4621 

(1995). 



